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ABSTRACT 



o 
Q 

The Cosmic Neutrinos Background (CNB) are Primordial Neutrinos decoupled when the Universe was 
very young. Its detection is complicated, especially if we take into account neutrino mass and a possible 
breaking of Lorentz Invariance at high energy, but has a fundamental relevance to study the Big-Bang. 
In this paper, we will see that a Lorentz Violation does not produce important modification, but the 
mass does. We will show how the neutrinos current velocity, with respect to comobile system to Universe 
expansion, is of the order of 1065 [^p] , much less than light velocity. Besides, we will see that the 
neutrinos distribution is complex due to Planetary motion. This prediction differs totally from the usual 
massless case, where we would get a correction similar to the Dipolar Moment of the CMB. 



> INTRODUCTION 

O 

From the beginning, the photons and all particles were coupled forming a plasma that was evolving 
under the influence of the Universe expansion. In such a moment, when the photons were dominating 
the expansion, the neutrinos were decoupling from the plasma and evolved in an independent way. One 
of the last discoveries about neutrino is its mass. This has relevant effects in the Standard Model and 

j. in some of its characteristics, distinguishing it from the photons. One of them, which we will study, is 

its velocity. Thus, we will analyze the evolution of the neutrinos's kinetic energy since its decoupling till 
today. 

Other phenomenon that we will study, that is directly related with the first one, is the neutrinos 
distribution. The detection of the Cosmic Microwave Background of Photons (CMB) is the best proof of 
the Big-Bang scenario [T] that helped to check or refute models that describe it, and study the composi- 
tion of the Universe. Because of this, it is important to study the Cosmic Neutrinos Background (CNB), 
especially the form of their Distribution Function to consider the effect of the peculiar velocity of the 
planet, named Dipolar Moment in the CMB, and optimize the detection. This is already complicated 
due to the low interaction that the neutrinos have with ordinary matter. The calculation will be done 
for photons and neutrinos in parallel. 

Finally, we will include a Lorentz Invariance Violation (LIV) represented by an alteration to the 
Dispersion Relation of energy given by [21 El SI [5] : 
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ui2 2 2,2 4 

E = V maxP + ™ C 

Where v max = c(l — a) is the maximum attainable particle velocity with a ~ (10~ 22 — 10~ 23 ) . The 
motivation to use this LIV comes from the possibility that, at the high energies available in the Big Bang 
there take place some LIV due to Quantum Gravity [2111115]. If such a LIV exists, the first problem is the 
appearance of a privileged reference system, but fortunately exists a natural candidate, the one where 
the CMB is isotropic. A LIV without a preferred frame as in Double Special Relativity [5J, will not be 
considered here. 



1 ENERGY AND VELOCITY OF PRIMORDIAL NEUTRI- 
NOS. 

Initially the neutrinos were in thermal equilibrium with the rest of matter. For this, is necessary that 
Ti 3> H, where L; is the rate of interactions of the species i, H cx T 2 is Hubble's constant and T the 
temperature. While the neutrinos are kept in equilibrium, its distribution will be given by Fermi-Dirac's 
statistics: 

f eq (E,T) = 

e k B T + l 

During cosmic expansion, the temperature will be diminishing down to a point where T u < H and 
Yi^v 3> H. This means that the neutrinos lost the equilibrium and are decoupled from the rest of the 
matter. We will name T Ut o the neutrinos decoupling temperature that are obtained when we impose 
Y v ~ H(T). To see what is happening with its distribution we will do the following analysis. For a time 
to, an observer sees in any direction a quantity dN = fdrdp of neutrinos in an volume d 3 r and with 
momentum between p and p + dp. After a dt time, the neutrinos have not interacted, so dN remains 

constant, but the volume in which they are, have increased in a factor (^^j§j^y^j an d the momentum has 

diminished in gz|^p^ , because of the expansion of the Universe. This means that f(E,T v ) is constant 
in time. Therefore, for t > to (or T„ < X^d) with to the moment in which is produced the decoupling, 
the distribution function is given by [7]: 



f[E(p(t)),T v (t)] = f eg [E(p D ),T u , D ] = /, 



eq 



(i) 



When the subscript D refers to the age of decoupling. In addition, we know that the number of 
neutrinos, the total energy and the energy per neutrino are given by: 

Nu = -*rLi I f(p,T u )d 3 P (2) 



(2-Kh) 

E »= 7^ I E(p)f(p,T„)d 3 p (.'-!) 



(2irh) 



Where: 

E 2 (p)=v 2 maxP 2 + m 2 c 4 
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to allow for a small LIV in the dispersion relation. 



Now we can determine the Distribution Function that they will have after being decoupled. It is 
possible to express the energy of the neutrinos (high energies and small masses) as E(t) = v max ^p{t) 
during the decoupling (We use an expansion with zero order in the mass because / depends exponentially 
on E), and as pu = P(t)^~ we obtain using (jlj: 

f\p,T v ] = -^z (5) 

e fc s T - + 1 

With T v = T„ jjt^jj and /i„ — because of the low interaction that they have with matter. This means 
that the distribution of neutrinos after decoupling is Fermi's with temperature T„, therefore RT V = etc. 
Replacing it in ^ and ^ : 



(2tt/i) 3 J e kBTv + j 

Naturally, N v will be constant in time. Using the change of variable x = "^y 81 , we obtain: 

N » = ( 6 ) 

Where £(3) = 1.2021 is the Riemann's Zeta function. We see that, in fact, N v keeps constant in 
time because V oc R 3 (t) and T v oc i? _1 (i). To determine E v , we must compute the integral, which is 
complicated for the general case. Thus , we will analyze the extreme cases where the neutrinos continue 
being relativistic and when they do not. Due to the spherical symmetry, the velocity is only radial, 
therefore we just must determine its modulus. The modulus of the velocity of a particle is given by: 

de 
dp 

Being e and p the energy and the momentum of a particle, related by our dispersion relation: 

2 2 2,2 4 

£ = v rnaxP +m C 

While the particle continues being relativistic, developing the derivative till the second order in the 
mass (e 3> mc 2 ), we obtain: 

/ 1 /mc 2 \ 2 \ 
v p ^ v max ,u I 1 ~ 2 V ~e~ J I 

Notice that we must use E(p) — v max p to calculate E v , to the order of approximation in the mass 
that we are considering. 

Now, if the particle becomes Non-Relativistic, we have that the energy and the velocity of a particle 
to second order in the momentum (pv maXil , <C mc 2 ) will be: 

2 , ( Vmax,v \ 2 P 

£,, ~ mc 1 ' 1 



2m 



^-(^) 2 ^-W,J 2 (^-l) (8) 
op V c / m V \mc z / 
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Where we see that, to keep the order in the momentum, the calculation must be up to second order 
in the expression of E(p), therefore E{p) — mc 2 + C""'" ) 



1.1 Relativistic Neutrinos 

As we said, to determine E v we must use E{p) — v max p. With this, we obtain the expression: 

" 1'10/rv.':. 



Using (|6| and (JoJ in @ and 0, we obtain: 

77T 4 



180C(3) 



k B T v (10) 



1 /180C(3)m„c 2x 2 " 



v v = v max , v \l-^ j j (11) 

If we define the relative velocity between the neutrinos and the photons as Av — c — v v , result: 

, v max v ( 180((3)m u c 2 



„2N 2 



Where Av max = c — v maXlV = ca v . We can see that this factor vanishes if the violation does not exist. 
Evaluating numerically: 

Where we have separated the LIV dependent part from the rest. 
1.2 Non-Relativistic Neutrinos 

In this case we have that E{p) = m v <? + ( Vma *-" ) 2 ^L, therefore, when we evaluate in E v using (6), we 
obtain: 



*-w(i + i(ig)'£) <i3) 

With I n = J °° -^f^dx = (l — ~) n!£(n + 1). Then, evaluating in O and (Isl) , we have: 



C(5) fk B T v 



2\ 



LC(5) k B T v 

v » = v — V 15 c(3) MT (15) 



Giving a relative velocity: 
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- - ^ilM k -^ +(l- ^ k -ff] (16) 

c w y c(3) m v \ v c(3) m„ y ^ j 

Where we have separated the LIV part from the rest and £(5) = 1.0369. Since the neutrino velocity 
cannot be higher than its maximum velocity, we must see to what temperatures this approximation is 

valid. We have that v v > v maxM if kgT u > ^?^ M V . It means that the approximation is valid if 
k eT v « \j$fe)M v ~ 0.28M„. 



1.3 Numerical Results and Analysis 

In the age of decoupling of the neutrinos, we know that ksT V £, ~ (2 — 4) [MeV] and, currently, /cbT^o 
1.68 x 10 -4 [eV]. In addition to this, for cosmological parameters, we know [5]: 



m v . < 0.17[eV] 



That clearly indicates that they are relativistic in the moment of the decoupling. There exist many 
estimations of the masses of the neutrinos, but none of them are very precise. Thus, we will use 
m v ~ 0.17[eV]. This way, we are sure of being inside the correct limits and we will find the maxi- 
mum effect that the mass could have in the velocity of neutrinos. This way, none of these estimations is 
below ksTv o, therefore they are Non-Relativistic nowadays. 

Before discussing the results, we will analyze the effect of the LIV. Thus, we will compare our rela- 
tivistic expressions with our Non-Relativistic ones. If we observe these expressions, we can see that both 
are proportional to v maXjl/ , which is the only thing that depends on a v . It means that the difference in 
percentage between the case with and without LIV is always: 

^(°)-^K) 100% = avlm% = 1 x 10 -20 % 

Therefore, it is not possible that this LIV has an important effect in the neutrinos, then we will 
continue our calculations using a = 0. 

Previously we mentioned that our Non-Relativistic approximation is valid if k fJ v < 0.28. At present 
we have that k f^ v — 10~ 3 fulfilling the Non-Relativistic bound, but with a mass 100 times minor 
(~ 2 x 10 -3 [eV]) the bound is not respected. However it does not correspond to the relativistic case 
either. To be kept relativistic, we need a mass 10000 times smaller or less (~2x 10~ 5 [eV]). 

In Figure [3] the evolution of the velocity of the neutrinos due to the expansion of the Universe is 
represented graphically . We define the adimensional quantities z = fe ^, , y = It is indicated in the 
graph that the time grows towards bigger values in z. Clearly, we see that the neutrinos suffer a rapid 
deceleration from the time of decoupling . Then, this deceleration begins to diminish slowly, approaching 
a zero velocity. 

All the estimations of M v indicate that we are in a zone dominated by the Non-Relativistic approxi- 
mation. The estimation for the smallest masses ranges between 10~ 4 and 10 -3 [eV]. Remembering that 
our top limit is 0.17 [eV], we see that we are currently in the region 0.6 < z < 1012, which is a very wide 



range. Evaluating numerically in (15 1, we obtain v v = 3.55 x 10 3 c = 1065 [— ], with a mass of 0.17 
[eV]. This velocity will be bigger if we use smaller neutrino masses. 
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Up to now, we have assumed that the neutrinos are not affected by the galactic potential, they are 
free particles and are not relics from the Milky Way [HJ-PZI]- To check this point, we consider the relation 
between kinetic and potential energy of the neutrino in the Milky Way. That is: 

m v v 2 GMm v 
2 = R 

/2GM 

Where v would be the limit velocity where the potential energy is comparable with the kinetic energy. 
Evaluating in M ~ 2 x 10 42 [kg] and R ~ 4.7 x 10 20 [m], mass and radius of the Milky Way respectively, 
we obtain v ~ 754 [— 1. Since v v > 1000 [— 1, our supposition is correct. 



2 THE CNB DISTRIBUTION. 

To determine the effective neutrinos distribution (distribution from Earth), we need to use equation Q 
in the comobile system of the Universe. In addition to this, we will use the photons distribution when 
they are decoupled, that is: 

f(p,T) = ^J— 

ekT — l 

That corresponds to an ultra-relativistic Bose-Einstcin's distribution with RT = cte. Currently, it 
has a temperature of 2.73 [K]. In addition, we saw in the previous chapter that a LIV of the form: 

ui2 2 2,2 4 

E = v max p +m c 

Is not markedly different in its energy and velocity in comparison to the usual Dispersion Relation 
(vmax = c). Because of this the usual special relativity rules are valid. For instance, Lorentz's Trans- 
formations for the coordinates of space-time and of energy - momentum. Then, we can compare the 
neutrinos distribution in the comobile system and Earth. It is possible to demonstrate (See Appendix): 

f'(p',T>)=f(p,T) (17) 



E = j(E' - v t p' cos (00) (18) 

Where the primed elements refer to the reference system of Earth and the non primed to the comobile. 
& is the angle that is formed between the vision line and the direction of Earth motion and Vt is the 
Earth's velocity |llj . We can see that the distribution function is invariant under Lorentz's Transforma- 
tion and the energy changes with the angle of vision. 



Now we will analyze some cases. First, the photons of the CMB to guide us because they are already 
very well known , and secondly, the neutrinos. Using expression (18) we will determine p' as a function 
of p. 



G 



2.1 Photons 

In this case, we have that E = cp, therefore the expression ( 18 ) is reduced to: 



l-^cos(0') , 
P= I 



Replacing in (17 1, we obtain: 



/ (f> , r 7 ) = / | -p , T 7 



As the photons, after being decoupled, continue with a distribution of the form: 

A 1 



We can leave our expression as: 



/'(p',T 7 ) = / |p',T 7 



1 - 2* COs(0') 



Therefore, the photons distribution detected from Earth, /', in a specific direction, will be of the same 
form that the one detected in the comobile system to the Universe, but with a different temperature given 
by: 



rpl rp \l ( c ) 



7 7 1 - cos(e') 

If we consider that «t <C c, we have: 

T 7 ~ T 7 (l + ^ cos(0')) 



Ar 7 _ «t 

T 7 c 



cos(0') (19) 



that is known as the Dipolar Moment, and is of the order of 10" 



-4 



2.2 Neutrinos 

Now, we have particles with mass. Currently, the neutrinos ar e N on-Relativistic, therefore E = m u c + 

2 I I 

2^- for both the comobile and Earth systems. Evaluating in (18) and using the approximation v t <C c 
up to second order in p' and Vt, we obtain: 

p 2 = p' 2 - 2m v v t p' cos(8') + mlv 2 



Evaluating in (17 1, we have: 

f'ip', K) = f Up>*-2m v v t p'cos(6>)+mlvl T v 



(20) 
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a(8>) 


Direction of Observation 


1 


In favour of the Terrestrial Movement 


0.5 


60° deflected to the Terrestrial Movement 





Perpendicular to the Terrestrial Movement 


-0.5 


120° deflected to the Terrestrial Movement 


-1 


Against the Terrestrial Movement 



Table 1: Directions of Observation. 

Values of a{9') used in Figure R] with the corresponding direction of observation. 



In this case is impossible to find a relation between T' v and T v , but we know that the distribution is 
given by Seemingly, we can only notice the effects graphically. To facilitate our analysis, it will be 
helpful to define the number of neutrinos per solid angle d£l' of momentum as: 

dN gV , , , , 2 
^ = (2W /(2, ' T> dp 

With this, we can obtain the distribution function of the number of particles: 



Fyx) = j£^ap'X)^ (2D 

In our case, the distribution function F 1 will be: 



„/2 

F'ip'X) oc P 



+ 1 

2.3 ANALYSIS 

To do our analysis, it is useful to introduce the adimensional variables x = t ^ , a(8') = ™' B ^ C cos(9') 
and b = a(0). With these parameters, our distribution is: 

9 

F' <x ; X (22) 

Considering a terrestrial velocity v t ~ 300 [— ], we can see that b ~ 1 for M v — 0.17. It means 
that —1 < a(8') < 1. This range will be smaller if we use a smaller mass, but then the Non-Relativistic 
approximation is less precise. 

In Figure [3] it is shown F'; here we have used a value of b ~ 1 and some representative values of a(9') 
(See Table 0. 

Let's remind that the distribution F' represents the particles number that come from certain direc- 
tion and momentum. We see in Figure [4] that the distribution suffers a loss of homogeneity, which is 
translated in more neutrinos observed in favour of the Earth's movement, but simultaneously the form of 
the distribution function is altered much more with regard to the distribution of the comobile system. If 
we move away from this direction, the neutrino number detected diminishes considerably and the small 
momentum are favored. 

The distribution maximum must fulfill the equation: 
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{^\/~^rnax 2aX maa; -j- -^max ^mfllj g\/ ma;c maT+ _j_ x^- iax 2ttX maa: -j- 6^ — (23) 

It is complicated to find a general expression for xm<ix, although a numerical treatment is readily 
available. As an example, we can study the extreme cases a{9) — b and a(6) — —b where b ~ 1 for 
M v = 0.17 [eV]. Evaluating in we obtain: 

Xmax(a = b) = 2.463 
Xmax(a = —b) = 2.091 
This means that the momentum of the majority of detected neutrinos will be: 

2.091 < < 2.463 

This will be a useful information to plan the detectors. 



Now, if we use a smaller mass, the differences between different directions in the distributions di- 
minish. In Figures [5] and [6] we can see the distributions with a mass 10 and 100 times smaller, where 
the Non-Relativistic approximation can be still valid. To compare, the photons distribution appears in 
Figure [7] for different observation angles. Comparing this with Figure |6j we see that the effect produced 
in the neutrinos distribution is bigger always than the produced in the photons. 

This happens because the photons always go to a bigger velocity than the terrestrial (c^> v t ), there- 
fore Earth would be almost still with regard to the comobile system, doing that the isotropy almost does 
not change. On the other hand, if the neutrinos acquire mass, they will be found submitted to a decel- 
eration as the Universe expands (See Figure |3| so that that currently the neutrinos are Non-Relativistic, 
with a velocity not much higher than vt ■ This means that the effect of the terrestrial movement begins 
being important in the velocities addition and will be increased with the time due to the constant cooling 
of the neutrinos. This will reach the point in which the neutrino velocity will be much smaller than vt 
and, practically, the planetary movement will predominate. This will be reflected in an increase of the 
distribution in the direction of the terrestrial movement. 



CONCLUSION. 

The mass of the neutrinos brought important modifications to its velocity. Without mass, the neutrinos 
would have supported a constant velocity and equal to the light velocity, c. On the other hand, with 
non zero masses, its velocity is affected by a strong deceleration (See Figures |3|, therefore they are Non- 
Relativistic nowadays. As we have developed an expression for the velocity with regard to the comobile 
system to the expansion of the Universe, it is necessary to use the addition of velocities to determine the 
mean neutrino velocity relative to Earth. Thus, we use Lorentz's Transformations since the LIV did not 
bring any important effect. The difference that is produced in its velocity with and without LIV is of 
~ 10~ 20 %, which is totally negligible. Then, we can use the invariance of the distribution function to 
relate the comobile system to the terrestrial. 

In the same way, the mass of the neutrinos brought important changes to the distribution. Unlike the 
photons, it was not possible to introduce a similar term to the Dipolar Moment because the temperature 
would depend on p' . Greater the mass greater the effect. In addition, the distribution is widely favored 
in the Earth's direction, but if we move away from this direction, the neutrinos number diminishes. In 
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spite of that the variation depends greatly on the mass; as time goes by, the neutrinos will be cooling 
diminishing little by little its velocity This means that in some moment the velocity of the neutrinos 
will be less than the terrestrial speed. In the future, the neutrinos will be almost still in comparison to 
the Earth's velocity. In this moment, we will only detect the neutrino that "crash" with Earth when it 
advances. 

To sum up, we see that the existence of the neutrino mass produces a relatively important effect in its 
evolution, which is reflected in the perception that we have of them especially in the loss of homogeneity 
in the distribution function. Thus, for its detection is advisable to use detectors of neutrinos directed in 
favour to the terrestrial movement or to use a satellite located in someone of Lagrange's points of the 
Solar System to keep the isotropic distribution of the comobile system, as the satellite Planck Surveyor 
that will observe CMB [12] . 
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APPENDIX. 

The special relativity rules say us that two reference systems can be related by: 

x = x'j_ + 7 (at' 1 1 + v t t') 



and 

P = Px + 7 \ + 
E = j(E' + v t ■ p') 

Where the primed reference system is moving away from the non primed to a velocity v. The coeffi- 
cients with the subscripts |j and _L represent the parallel and perpendicular components of the velocity 
v respectively. Since we are considering a particle in the universe, our primed and non primed reference 
systems will be, respectively, Earth and comobile system to the Universe expansion, therefore v is the 
planet velocity. Thus, from now, we will call its v*t- 

If particles go to Earth along the vision line (See Figure [TJ, the Lorentz's transformation can be 
written as: 

( v tx\ \ 

x x ,i = x' ± j x l{ = 7 (4 + v t t') t = 7 -1 (24) 



P±.i=P±.i P\ t 

c 



-y(p'\\-^ E ') E = 1 {E'-v t p\) (25) 



Where i = 1, 2 label both perpendicular coordinates to Vt- Its differential form considering an 
instantaneous measurement from Earth, is t' = cte or dt 1 = 0, giving: 

Vtdx'u 

dx± i = dx\ i dx\\ = *ydx'u dt = 7 — = — (26) 

1 11 c i 

dp ±<i = dp' ±4 dp\\ = 7 (dp\ - ~dE'^ dE = j(dE' - v t dp\ ) (27) 
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s s" 




Figure 1: Description of both reference systems. S: Comobile reference system to the Universe expansion. Earth has a 
velocity vt and the neutrino has momentum p. Between both we have the vision angle 6. S': Earth reference system. Earth 
is still and the neutrino has momentum p' . We have the vision angle 9' measured from Earth. The coordinates system of 
S and S' are related by the Lorentz's Transformation. 



Now, when we count the particles number from Earth in a specific direction instantaneously (dt' = 0), 
inside a volume d 3 r' we have dN particles with momentum between p' and p' + dp'. In addition, we know 
that dN is given by: 



dN = f'{p',T')d 3 p'd 3 r' 



(28) 



Where /' is the distribution function on Earth. In the comobile system, the particles are in a volume 
d 3 r and with values of momentum between p and p + dp, but in a dt time, given by (26 1 (different from 
zero because dt' = 0) some particles enter or exit of d 3 r. Thus, the particle number, in this system, is 
given by: 



dN = f(p, T)d A pd A r + f(p, T)d 3 pdS ■ udt 



(29) 



Where u — c 2 ^ is the particle's velocity and dS is the differential area, with normal direction. Both 
expressions for dN are, simply, a variation of the continuity equation: 



dp 
di 



+ V • (pu) = 



With p = f(p,T)d 3 p. Since dN must be the same in both systems, we must equal (28 1 and (29). 
Then: 



With (See Figure [2]): 



dS 



f(p>, T')d 3 p'd?r' = f(p, T)d 3 p d 3 r + c 2 dt^- 



d 3 r = dxu A dx_\_ t \ A cfe_L,2 
-{dxn A dx±,i£± t 2 + dxu A dxx^&x^ + dx± t i A dx± i 2X«) 

p= -(P\\X\\ + P±,lX± t l + P±,2X± i2 ) 



(30) 



Where A represents the anti-commutative product between the differentials. Evaluating in (30 1, using 



(251 and (26), we have: 
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X. 2 



Figure 2: Representation of the volume element d 3 r, where we see the surface elements. We can see, clearly, the vectorial 
direction of p and dS. 



/V.^WdV = f(p,T)d 6 pd d r 



~E 



(31) 



Replacing (27 1 in d 3 p = dp± t i A dp± t 2 A dp», we obtain: 



d A p = dp' x>1 A dp' ± 2 A 7 - ^|dE' 

But we know that £" 2 = c 2 (p'^ + 2 + f>'||) + m 2 c 4 . Deriving, we obtain the relation E'dE' 
c 2 (p'j_ idp' x 1 + p'j_ 2 dp± 2 +P\\dp\\)- Evaluating: 

d 3 p = dp' ±>1 A dp' ± 2 A dp'ff 1 1 - j 



Where we have used that dp' ± , A dp' ± { = for anti-conmutativity. Using (25), d 3 p stays 



E 



d3 P = d3 p'^, 



Then, (31 1 is reduced to: 



f(p',T')=f(p,T) (32) 

This means that the distribution function is Lorentz invariant. In reference [T3], this has been discussed 
differently. They used: 



dN = f(p,T)d 3 pd 3 r 
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Naturally, they obtained that / is not Lorentz invariant. 
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Figure 4: Primordial Neutrinos Distribution in the current age (Ec |22| for different values of a and b ~ 1 that corresponds 
to M v = 0.17 [eV]. The black curve represents to the distribution in the comobile system. 
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Figure 5: Primordial Neutrinos Distribution in the current age (Ec |22| for different values of a and b ~ 0.1 that corresponds 
to M v = 0.017 [eV] 
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